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SUMMARY 


An investigation was made of a thin-walled cylinder under axial 
compression and various internal pressures to study the effect of the 
internal pressure on the compressive buckling stress of the cylinder. 

A theoretical analysis based on a large -deflection theory was also made. 
The theoretically predicted increase of compressive buckling stress 
due to internal pressure agrees fairly well with the experimental 
results. 


INTRODUCTION 


The buckling of thin-walled cylinders under axial compression and 
lateral pressure has been investigated by Flugge (reference l) who 
found that the effect of the internal pressure on the buckling load is 
negligible. Flugge' s conclusion is in contradiction to the results of 
a series of tests, made at the .Langley Aeronautical Laboratory of the' 
NACA, of two curved panels under axial compression and various lateral 
pressures. These test results, reported in reference 2, showed an 
appreciable strengthening effect of the lateral pressure on the buckling 
load of the curved panels. The apparent discrepancy between these 
experimental results and the prediction by Flugge .'s theory made it 
desirable to investigate this problem further. Consequently, additional 
tests were made of a cylinder under axial compression and various 
internal pressures for which results are presented herein. A theoretical 
analysis of this problem is also presented which differs from that of 
Flu'gge in that the present analysis is based on large, rather than small, 
deflection theory. 



2 


NA.CA TN 2021 


APPARATUS AND PROCEDURES 


Test specimen .- Thei specimen used for the tests was a cylinder/ 

32 inches long with a 15 -inch inside radius, made of 24S-T aluminum 
alloy sheet of 0.0249-inch average thickness. It was closely riveted 
around two heavy steel rings, one at each end. The butt-joint of the 
two longitudinal edges was covered both inside and outside by straps' 

0.032 inch thick and 1^ inch wide, along the total length of the cylinder. 
(See fig. 1. ) 


The two heavy steel rings were made of — - by 4-inch steel bar stock 
rolled to the diameter of the cylinder. Two hy 2-inch spreader bars 
were used to reinforce the ring as shown in figure 1. A ring with a 
flange, machined flat, was fastened to the by 4-inch steel ring to 


provide an even bearing surface on which a steel cover plate was fitted. 
Three steel blocks were placed on top of the plate. The applied 
compressive load was transmitted from the machine head through the three 
steel blocks to the cover plate. The joint between the cylinder and the 
cover plates was sealed. 


Equally spaced along the inside circumference of the cylinder at 
midlength were 16 strain gages, and directly opposite to them on the 
outside were 1 6 more gages. These gages were placed to measure strains 
in the longitudinal direction. Six more gages, three inside and three 
outside, were placed to measure the circumferential strains. 


Test procedures . - The specimen was subjected to compressive load 
in the 1,200, 000 -pound universal testing machine of the Langley 
structures research laboratory. Compressed air was used to produce 
internal pressure, which could be maintained at any desired constant 
value. The pressure was measured by a manometer. The strains were 
recorded by standard electric strain -gage equipment and the end- 
shortening was measured by dial gages. 

The cylinder was preloaded and the strain-gage readings were 
taken. The three steel blocks were so adjusted that all longitudinal 
strain-gage readings around the circumference of the cylinder were equal. 

The compressed air was then let into the cylinder unt il the 
desired internal pressure was reached. The axial compressive load 
was increased in increments until buckling was observed. At each load 
increment, all gage readings were recorded. The load was then decreased 
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until the buckle a disappeared and increased a second time to check the 
reading obtained the first time. During all these steps the internal 
pressure was maintained constant. 

The axial load was then reduced and the internal pressure was 
changed to another value. For each value of internal pressure the same 
procedure was repeated. 


EXPERIMENTAL RESULTS 


A typical experimental result is shown in figures 2(a) and 2(b) 

for the case in which the internal pressure was 1^ psi. In figure 2(a) 

the compressive load is plotted against the strain-gage readings for 
four different pairs of gages, within the range where the load-strain 
relatioii is linear. In figure 2(b) the strain-gage reading is plotted 
for all strain gages at three compressive loadings close to the buckling 
load. Figure 2(b) indicates that buckling occurs at a compressive load 
of 12,700 pounds between strain gages 22 and 23. (Note the intersection 
of the curves at two consecutive loadings.) A buckle at this location 
was observed during the test. The compressive load at which this 
phenomenon occurs is considered the buckling load. 

Since the buckling occurs locally and not simultaneously at all 
the gages, the local buckling strain is obtained by dividing the buckling 
load by the slope of the linear portion of the load-strain curve corre- 
sponding to the gage at which the buckling occurs. The corresponding 
stress is the buckling stress. The buckling stresses for various 
internal pressures were determined in this same way. 

The results are tabulated in table 1 and plotted in figure 3 In 


the internal pressure, R is 

thickness, and E is 
Young's modulus. Except for the first test corresponding to p = 0.1028 
in which the cylinder had undergone no previous buckling, all the tests 
were carried out on the cylinder with possible permanent set. ' 


terms of the two nondimens ional parameters 


J cr “ E 


u cr R 
t 


- e/r 

P = flt 


where a Ucr is the buckling stress, p is 
the radius of the cylinder, t is the wall 
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THEORETICAL RESULTS 


A theoretical analysis for calculating the buckling stress of a 
cylindrical shell under axial compression and Inter nal pressure was 
obtained by a "large-deflection" theory for which details are given 
in the appendix. The large-deflection theory was first advanced by 
Yon Karman and Tsien (reference 3) in the study of buckling stress of 
cylindrical shells under axial compression (but without internal 
pressure). This theory was subsequently improved by Leggett and Jones 
(reference 4). In reference 3 the buckling stress was shown to depend 
on whether the load was applied by a rigid loading machine or by a 
dead -weight machine. In the present analysis, the loading machine is 
assumed to be rigid. 

The existing procedures for computation of the buckling stress by 
large -deflection theory involve the solution of four simultaneous 
nonlinear equations for each pressure loading. The numerical work 
is quite lengthy . The method used in the present study introduces a 
fifth equation which governs the conditions at which the buckling 
occurs. The fifth equation is based on consideration of conservation 
of energy, which is an extension of Tsien' s buckling criterion given 
in reference 5* Although a solution of five simultaneous equations 
is now necessary, the numerical work Is actually reduced to a small 
fraction of that required if the existing procedures were used. This 
reduction in labor is made possible through a proper choice of the param- 
eters, in the equations and the process of the computations. The results 
calculated by the present method are presented in table 1 and are 
represented by the solid -line curve in figure 3- The curve is cut off 
at a value of 5 cr = 0.605, corresponding to p «5 O.169. This constant 
value of 3 cr = 0.605 for p > O.169 is the same that is obtained by 
the classical' theory . 


DISCUSSION AND CONCLUSIONS 


From the theoretical and experimental results shown in figure d, 
the internal pressure is seen to have an appreciable strengthening effect 
on the cylinder. Although the two curves obtained from theoretical 
and experimental results do not coincide, both. show the same trends 
as regards the effect of internal pressure on the buckling stresses. 

If the increment of the buckling stress A3 cr due to the presence 
of internal pressure (that is, the difference between the buckling 
stress with the pressure 5 cr and that without the pressure (5 cr )- _ q) 
is plotted against the internal pressure, as shown in figure 4, a good 
agreement is -obtained between the theoretical and experimental results. 
These data indicate that, although further improvement of the theory 
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is necessary for the determination of the magnitude of the compressive 
"buckling stress, the theory gives a fairly good prediction of the 
increase of compressive "buckling stress that may "be expected as a result 
of internal pressure. The discrepancy between the theoretical curve 
and the experimental curve of figure 3 is believed to be caused by such 
factors as manufacturing imperfections in the specimens, material 
irregularities, and energy absorbed by the loading machine, which have 
not been included in the theory. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 

THEORETICAL ANALYSIS OF BUCKLING LOAD OF CYLINDRICAL SHELLS 
UNDER AXIAL COMPRESSION AND INTERNAL PRESSURE BY 
LARGE -DEFLECTION THEORY 
Background of Theory 


The use of large -deflection theory for shells under axial compression 
was first advanced By Yon Karman and Tsien (reference 3) in. an attempt 
to explain the discrepancies Between the Buckling- loads predicted By 
classical theory and those obtained from experimental results. (See 
for instance, reference 6.) The results of reference 3 indicated that 
cylindrical shells can Be maintained in equilibrium in the Buckled 
state By a compressive load considerably lower than that predicted By 
classical theory. A plausible explanation of this result is that, Before 
the classical buckling load is reached during a test, the cylindrical 
shell "jumps" from an equilibrium unbuckled state to an equilibrium 
Buckled state. The physical phenomena of the jump were further examined 
in reference 5 By Tsien. 

The treatment of Yon Karman and Tsien in reference 3 was left 
incomplete, however, in that the equilibrium positions at the Buckled 
state were determined By differentiating the total potential energy 
with respect to some But not all of the physical parameters involved. 

The resulting equations gave a relation Between the average compressive 
stress a and the end -shortening € in terms of the remaining parameters 
A set of curves of 0 against € were thus obtained for various 
combinations of the remaining parameters . 

Improvement of Yon Karman and Tsien ? s theory was made By Leggett 
and Jones (reference 4), who took the derivatives of the energy with 
respect to all the parameters and thus obtained a single curve 
Between a and €, representing all equilibrium positions of the 
cylindrical shell in the Buckled state. The same result was obtained 
By Michielsen (reference 7 ) in a similar process. Such a curve is shown 
By BC of figure 5 . 

Theoretically, when the cylinder is compressed, the -relation 
Between o and € follows the straight line ODA which represents 
the unbuckled state and will reach the point A if everything is perfect; 
the cylinder then Buckles and the relationship follows the curve ABC 
which represents the Buckled state. Before point A Is reached, however, 
some external disturbance may possibly cause the cylinder to jump from v 



MCA TW 2021 


7 


the unbuckled state represented by the point D to the buckled state 
represented by the point E. The positions of D and E on the respective 
curves depend on the actual physical conditions of the jump. 

' If the physical condition which governs the jump is known or 
defined, the buckling stress corresponding to the point D can be 
obtained directly without going through the labor of finding the curve 
ABC. This procedure can greatly reduce the amount of numerical work. 

In reference 5, Tsien introduced a criterion which governs the 
jump BE for the condition of loading obtained in a rigid testing machine, 
namely, that the strain energy, remains the same before and after the 
jump and that the jump occurs at constant end -shortening . According 
to this criterion the line DE must be vertical and must cut the curve 
AB in such a way that the two shaded areas ADG and GEE are equal. In 
fact, the area ADG represents the additional energy that is needed to 
assist the cylinder in jumping from a condition represented by D to that 
represented by G and the area GBE represents the energy that is given 
up by the cylinder when it arrives at- the lower energy level, point E. 

The energy represented by the area ADG is very small, and therefore 
a slight disturbance from the surrounding air might assist the cylinder 
to jump from the unbuckled state to the buckled state at a compressive 
stress well below the classical buckling stress corresponding to point A. 

Since the external disturbance is required to assist the cylinder 
to jump from the state corresponding to D to that corresponding to G, 
a slightly larger external disturbance can well cause the cylinder to 
make the transition from the state* represented by D* to that represented 
by B, except that in the case in which the cylinder jumps from D 1 to B 
the cylinder absorbs the energy of the external disturbance and does not 
reemit it. The buckling stress can be then as low as point D’. This 
fact was pointed out by Tsien in reference 8. 

In addition to the two criterions just mentioned, there are still 
others that might be used. In view of the fact that the choice of the 
buckling criterion is a much less important factor in the determination 
of buckling stress than are such other factors as, for example, the 
initial imperfections, Tsien's criterion of reference 5, as represented 
by the line DE, is as reasonable as any other, and the choice of this 
criterion greatly simplifies the numerical work. 

Tsien's criterion of reference 5 cannot be applied directly to the 
present analysis, however, because with the presence of the internal 
pressure the strain energy is no longer the same before and after the 
jump. In addition, the criterion is applied herein in quite a different 
manner from that of reference 5. In reference 5, a series of values of 
wave number . n and aspect ratio 3 were chosen and the criterion was 
applied to each pair of values of n and 3; the pair of values of n 
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and 3 which gave the minimum value of buckling load vaa considered to 
correspond to the buckling condition. In the present analysis, since the 
variation of a with e can be plotted only as a single curve, this 
criterion need be applied only once for each internal pressure. The 
results correspond to the minimum-potential-energy condition. 

In the derivation of the present analysis, the basic equations in 
reference 3 are first extended to include the effect of internal pressure 
Tsien's criterion governing the jpp for rigid machine loading . 
(reference 5) is modified, and the buckling stress is finally obtained. 


Symbols 


A list of symbols follows. Most of the symbols used in the present 
paper are the same as those in reference 3 ) exceptions are the use 
of p for Poisson's ratio, \ for wave length, and 3 for aspect ratio 
of the buckled waves. .< 


^•a 

^•b 

f 0> f l> f 2 
m 

n 


half wave length in longitudinal direction 

half wave length in circumferential direction 

parameters used in deflection function 

number of waves in longitudinal direction within 
length equal to circumference of cylinder 

number of waves in circumference 


p 


internal pressure 


t 




u 


w 


thickness of cylinder wall 

coordinates measured in longitudinal and 
circumferential directions, respectively 

component of displacement of a point on median 
surface of shell in x-direction ' 

component of displacement of a point on median 
surface of shell in radial direction 


a 


measure of average circumferential stress per 
wave length in longitudinal direction 


€ 


end-shortening of cylinder 
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<j 





P 


E 

♦ 


B i, B 2 , . . . Bg 

®1 ' > ' ) • • • ®6 * 


D l> D 2> • • • D 5 

< B >p - £ 


average compressive stress 
aspect ratio of ‘buckled waves 
Poisson's ratio 
total potential energy 
Young ' s modulus 
strain energy 

certain functions of 3 

certain functions of p and 3 
(D represents the functions Di, D2,... D^) 


(Dip - 0 |j 


(D represents the functions D-^, D2, . . . D5) 


R 


W, 


W,- 


W- 


Wt 


radius of cylinder 
elastic extensional energy 
bending energy 

work done by applied compressive load 
work done by internal pressure 


Nondimensional parameters: 


P 





■. a R 
a " E t 
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Subscripts : 


o 

u 

cr 


(Wi, W 2 ; W 3 , and are defined in the same 
manner.) 


pertaining to buckling condition 
unbuckled state just prior to buckling 
buckling condition 


Derivation of Basic Equations 

Three basic equations are derived in the following analysis to 
include the effect of internal pressure. They are the expression for 
the total potential energy, the expression for the strain energy, and 
the relation between the end -shortening and the average compressive 
stress. 


In order to calculate the total potential energy, the work done by 
the internal pressure should be included in addition to the energies Wj_ 
W 2 , and W 3 , which are given by equations (25), (26), and (27) of 

reference 3 as: 



A 2 B 2 B 4 P 4 C 2 


+ eV 

(1 + 93 2 ) 2 


(9 ■+ 6 2 ) 2 


b 4 h 2 

16(1 + B 2 ) 2 
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Wc 


in:tx a x b 6 ( ; 


rr^y(lf nl1 { f i s [V + 1* 2 ) 2 - d + 


(1 + P^fjfg + (1 + P k )f 2 2 


( 2 ) 


W- 


~EtX a X b 


2(! - ^ 2 )(f) 2 + n 2 |(|X + 3 2 )^! 2 + |fi f 2 + # > 2 S 


2 M§ f 0 + iW 


(3) 


where a, b, p, and v have "been changed to X a , X b , 3, and p, 
respectively, to agree with the notation of the present paper and where 


and 


A = ±fl 2 rfi - jifl + t 2 J 

B = ^f x 2n2 

0 - + f 2 J- ^ 

B + f 2 ) 

0 . Jfin 2 ^! + f 2 j 

B = n 2 (^fl + f 2 J 2 


The work done hy the internal pressure is, for a complete wave panel, 


P^a 



-k 


JO 


wp dx dy 

Jo 
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The nega.tj.ve sign is introduced because the radial deflection w is 
considered positive inward. 

If the same deflection function given in equation (l 6 ) of reference 3 
is used, that is. 


w 

R 



+ 




/ mx 
Icos — cos 


EL 

R 


+ ^ cos 


2 mx 

R 



cos 



+ 



2 mx 

R 


+ cos 



w 


the work done "by the Internal pressure "becomes 


Wp = -^pRX a l b (f 0 + 


(5) 


If the total potential energy 


0 = W x + W 2 - W3 - Wj 


is differentiated with respect to f Q and the derivative is set equal 
to zero, the following expression is obtained 


f ° + T - + f i f a + f 2 2 ) - ■« - 1 


(6) 


Substitute this expression into equations (l), ( 2 ), (3), and ( 5 ) for Wj_, 
W 2 , w 3 > and w p> respectively, and the following equations are obtained: 


W-, 




2 ♦ * (g 2 


E AEt 


(e 

n^f-, 4 + B 2 f L 3f 2 + B 3 f 1 2 f 2 2 + B^f^ 3 + ^f 2 ^ 

n2 ( 2B1 * + ^) f l 3 + ( 4B ^ + 3i) f l 2f 
( 2B 4 + ^) f l 2 + K f 2 + ^2 


( 7 ) 
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w 2 




- 'si ( s 5 f i 2 * w 2 ♦ 


v 


i. 


2 Etx a x b 


= 4 


2 (e) + ^(e ) 11202 ^! 2 + f l f 2 + f 2 2 j 


W 




£—.8 



&) — 1 - n 2( £&\l ( Ip • 2 . f f . f 2 

Et/ °ME \ Et \Et ]lm 4* 1 f l f 2 + f 2 


( 8 ) 


( 9 ) 


( 10 ) 


where 


= 64 


kj t 3 4 ^ II . 3 4 , 3 4 . P 4 

O T ». ^ + 


B 2 - 16 


8 




4 (1 + 32)2 d + 9 P 2 ) 2 . (9 + P 2)2 


A + §1 


B 3 "A 


L (1 + 3 2 ) 2 (l + 9P 2 ) 2 (9 + 3 2 ) 2 

4 


+ — 


n £ 

2 (1 + 3 2 ) 2 " (l + 9P 2 ) 2 " (9 + 3 2 ) 2 




i — L 


(1 + 3 2 ) 2 


b k = 


5 6(1 - n 2 ) 


1(1 + 3 2 ) 2 + £(1 + 3 4 ) 
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Equations ( 7 ) to (10) may be expressed in terms of the nondimens ional 


parameters o y 

P, p, T\, 

£ , and W as 


*1 = 

4^5 2 + 2p5p 

+ p 2 j+ + ti(P(- D3) + £ 2 D^ 

(11) 

w 2 = 

n 2 ^ 


(12) 

W 3 = 

85 2 + 8p5p 

+ 5 q^ 2 Dp 

( 13 ) 

V 

8pap + 8p 2 


(14) 


where 

D i - + p + p 2 ) 

I>2 = + B^p B3P 2 + + — B^p^ 

D3 = ^2B^ + + 2p) 

% = (2B4 .+ + |(p + p 2 ) 

+ Bg(p + p 2 ) 

The nondimens ional total -potential -energy parameter fi is 
0 = ¥ x + W 2 - W 3 - W p 

= - 4(5 2 + 2p5p + p 2 ) - (a - ipjn^Di + 
^ 2^2 " t 1^ I) 3 + ^ 2d 4 + 


( 15 ) 
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The nondimensional strain-energy parameter if is 


¥ = % + w 2 

= 4(5^ + 2p5p + p^) + ^2^4^ - 

^ 3 d 3 + ^ 2 D4 + n 2 ^ D 5 


The relation "between the end-shortening e and the average compressive 
stress can "be determined from equation (23) of reference 3 by integrating 
thus, 


04 a 



where a/E, as determined from equation (24) of reference 3, together 
with equation (6) herein is 


a pR 
E " Et 

Therefore, the relation between e and 5 in nondimens ional form 
..becomes 

r = £ 

t 

= 5 + pp + 


(17) 
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Equations (15), (l6), and (17) are for cylinders in the buckled state. 
For cylinders in the unbuckled state, the corresponding equations are 


0 U = -^(^ u 2 + 2pa u p + p 2 ) 

( 18 ) 

♦u = Ma u 2 + 2pa u p + p2) 

(19) 

iu = 5 U + W 

(20) 


Equilibrium Positions of Cylinders in Buckled State 


The equilibrium positions of cylinders in the buckled state can 
be obtained by differentiating the total potential energy of equation ( 15 ) 
with respect to each of the parameters q, £, p, and 0 and by setting 
the derivatives equal to zero. Four simultaneous, nonlinear equations 
are thus obtained: 


M.o. - 

di) 

8 . 0 . - 


= 0 = - 


( 5 " ’ 2 ( ti ^ 2d 2 + ■ 2 i 1 2d 5 

r; - ^pjnDi - 2(ni) 2 n 2 + i.5(t|C)d 3 - d 4 - 3 2 d 5 

(5 - -i-5jn(i>i)p - W) 2 0> 2 ) p + (nb(® 3 ) p - (d 4 ), 

5 ? 


(20 


’) 2 (d 5 ) p 


> (21) 


M _ 0 - - 

dp - ° - 


(° “ ■ppj T l( I) l)p - (ilC) 2 (l>2)p + (qU(B 3 ) )3 - 


(B4) p - 4 2 (D 5 ) (3 + -ppqD-L 1 


11 

P 
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where 

0>l) o = J3 2 (l + 2p) 

P 2 

(Dg)p = B 2 + 2B^p + + 2Bhg3 

(D 3 )p = 2 ( B4 + ®) 

Wp = |(1 + 2P) 

(D 5 ) p = B 6 (l + 2p) 

o 

(Di) p = P 2 (-* + P + P 2 ) = 2D X 

(Dg) p = + B 2 ’p + B 3 'p2 + B 4 'p3 + iB 4 'p 4 

(D 3 ) p . = 2B 4 '(1 + 2p) 

Wp “ 2 V 

(D 5 ) p = B 5 * + B 6 *(p +P 2 ) 

and 

P^ 17 p^ p^ p^ 

8+4 (1 + p 2 )3 + (1 + 9P 2 )3 + ^(9 + p2)3 

} p*+ p^ p^ 

- (1 + P 2 )3 + (1 + 9P 2 )3 + 9 (9 + p 2 ) 3 

Li P 4 + P 4 + 9 P 4 
2 (1 + P 2 )3 (1 + 9P 2 )3 (9 + p 2 )3 

BI4. ' = — 

2 (1 + P 2 ) 8 
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B 5 - = i— |p?(l + 3p2) 

6(1 - n 2 ) 2 


B 6' = 


6(1 - n 2 ) 


■P" 


Let 



The four simultaneous equations (21) become 


5' = (tiC) 2 (2D 2 ) - (n^)D 3 + 2 ti 2 D 5 (22a) 


O' = (tiU 2 (2D 2 ) - (hS)(1.5D 3 ) + D 4 + q2D 5 


a = 


(nC) 2 (D 2 ) - (nU(D 3 ) p + (D 4 ) p + t! 2 (d 5 ) ( 


% 

(Di) ( 


o' = 


(ti^) 2 (D^) ^ - (il£)( D 3 )p + ( D 4)p + ti 2 (D 5 ) 


!>1 


(22b) 

(22c) 

(22d) 


Theoretically these four simultaneous equations can be solved for q , 
£, p, and (3 in terms of 5 for a given pressure. If they are 
substituted into equation (17) , a relation between end-shortening ? 
and the compressive stress 5 is obtained which represents all 
equilibrium positions at the buckled state . In fact, this solution 
is essentially that obtained by Leggett and Jones (reference 4) and 
Michielsen (reference 7) for cylinders with axial compression but no 
internal pressure. 


Practically, however, the solution of the four simultaneous 
equations (22) requires a long and tedious numerical process. If only 
the buckling stress is required, calculation of only one point on 
the curve of a against £ rather than the whole curve is necessary. 
This solution can be obtained by the introduction of one more equation 
which governs the condition at buckling. 
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Buckling Criterion 

In reference 5, Tsien gave the following criterion which governs 
the condition at buckling: That the strain-energy of the buckled 

cylinder is the same as the strain energy of the unbuckled cylinder 
when the cylinder is tested in the rigid testing machine so that the end- 
shortening does not change during buckling. This criterion is apparently 
established from cons'iderations of conservation of energy. Although 
other physical criterions can be used (for instance, see reference 8), 
the criterion of reference 5 was chosen and extended to include the case 
for which. the internal pressure is present. The choice of this criterion 
simplifies the numerical work. 

When internal pressure is present in the cylinder, work is done by 
the pressure during buckling. The strain energy in the buckled state 
is no longer equal to that in the unbuckled state, but 

* = * u + (23) 


where 

or 



is the work done by the pressure during buckling, 


ZVWp =I 


EtX a ^b 


/b\ 2 


px a 

"lA-h 

w 

1 


p(w - w u )dx dy 


J 

0 J 

0 




(24) 


Equation (l4) can be rearranged as follows: 


\ = 8(p 2 + nap) - -^pq^D-L (25) 

P P 2 

Therefore, for the unbuckled state, the last term is eliminated and 


(W p ) u = 8 (P 2 + P 5 uP) 


(26) 


Then, from equations (24), (25), and (26) 


AWp = 8np(5 


5 J 


■ i p71 ^ Di 

p 2 


(27) 
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The buckling criterion becomes (equations (23) and (27)) 


$ = + " 5 u) "‘^ 1 l^ 2l) l 

or, from equations (l6) and (19), 

U(5 2 + 2\io§ + p 2 ) + vpt,^T>2 - + ^ 2d 4 + n 2 ^ 2 I>5 

= 4(5 U 2 + 2|aa u p + p 2 ) + 8(5 - 5 u )|ip - -^PuS 2 ^ (28) 

* ^ 

Since the end -shortening remains unchanged during buckling, that is, 

€ = the following relation is obtained from equations (17) and. (20). 


5 u = 5 + iri^D-L 


(29) 


If this relation is substituted in equation (28) and if the relation 
= ~ use< t> the buckling criterion becomes 


5' = (ti0 2 (d 2 - Jg>! 2 ) ~ (t|C)D 3 + \ + q 2 D^ 


(30) 


The solution of the five equations (22a), (22b), (22c), (22d), and (30) 
gives the buckling stress for a given internal pressure. The following 
section presents a very simple method for the solution of these five 
simultaneous equations. 


Method of Solution 

From equations (22b) and (30) and equations (22a) and (22b ) , the 
following equations are obtained: 




»3 

2 (D 2 + ^D x 2 ) 

D4 - |d 3 (^) 
D 5 


(31) 


(32) 
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For a preasaigned value of P, assume various values of p and 
compute tj t, and q 2 from equations (31) and (32). Substitute these 
values in equations (22a) and (22c) to obtain (5') and (a') c ; 

respectively. Plot both (d') a and (d') c against p. The intersection 

of these two curves determines a pair of values o' and p which are 
called 5’ 0 and p Q . The corresponding (q£ ) Q and (q 2 ) Q are computed 

and substituted in equation (22d) from which the pressure p can be 
calculated. For each assigned value of P, there are obtained 
corresponding values of 3 ' Q and p. A curve of 5' against p 

can thus be determined. If the following relations are used, 

5 ‘o = (T D l) 0 ^o - 

^cr = (5 u )o ~ + g(Tj^Di) o 

the relation between a cr and p is obtained as shown in figure 3. 

Cut-Off Buckling Stress 

When equation (31) is derived from equations (22b) and (30), a 
factor (tj£) =0 is also obtained. If this relation is used instead 
of equation (3.1), it can be shown that the buckling stress a cr 

can never exceed the classical buckling stress 0.605 which is independent 
of pressure. 
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TABLE 1 

BUCKLING STRESSES FOB VABIOUS INTERNAL PRESSUBES 


Experimental 

Theoretical 

P 

a 

Aa cr 

P 

a 

A *cr 

0 

.01715 

.03425 

.0514 

.0685 

.0856 

.1028 

0.1936 

.252 

.277 

.309 

.350 

.363 

.407 

0 

.058 
.083 . 
.116 
.156 
.170 
.213 

0 

.02 

.04 

.06 

.03 

.10 

.12 

.14 

0.376 

.444 

.480 

.506 

.528 

.547 

.565 

.581 

0 

.068 

.104 

.130 

.152 

.171 

.189 

.205 
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Figure 3 *- Comparison of theoretical and experimental results of the 
buckling stress at various internal pressures. 



MCA TN 2021 


Theoretical results 
Experimental result, first test 
Experimental results, subsequent tests 





